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Three red giants with substellar-mass companions! 

A. Niedzielski, 1 A. Wolszczan, 2,3 G. Nowak, 4,5,1 M. Adamow, 6,1 K. Kowalik, 7,1 G. 
Maciejewski, 1 B. Deka-Szymankiewicz, 1 M. Adamczyk 1 


ABSTRACT 


We present three giant stars from the ongoing Penn State-Torun Planet Search 
with the Hobby-Eberly Telescope, which exhibit radial velocity variations that 
point to a presence of planetary -mass companions around them. BD+49 828 
is a M = 1.52 ± 0.22 M 0 K0 giant with a rn sin i = ITIq^Mj minimum mass 
companion in a = 4.2^;2 2 AU (259CU 3 4 5 6 7 gQd), e = 0.35+q; 2 q orbit. HD 95127, a 
log L/L 0 = 2.28 ± 0.38, R = 20 ± 9R 0 , M = 1.20 ± 0.22 M 0 K0 giant has 
a msmz=5.01^oi4 Mjminimum mass companion in a = 1.28^°;^ AU (482^d), 
e = O.ll^QQg orbit. Finally, HD 216536, is a M = 1.36 ± 0.38 M 0 K0 giant 
with a rn sin i = I^Iq^Mj minimum mass companion in a = 0.609^qqq 2 AU 
(148.6^'^), e = 0.38 ^q!q orbit. Both, HD 95127 b and HD 216536 b in their 
compact orbits, are very close to the engulfmcnt zone and hence prone to ingestion 
in the near future. BD+49 828 b is among the longest period planets detected 
with the radial velocity technique until now and it will remain unaffected by 
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stellar evolution up to a very late stage of its host. We discuss general properties 
of planetary systems around evolved stars and planet survivability using existing 
data on exoplanets in more detail. 

Subject headings: planetary systems-stars: individual (BD+49 828 b; HD 95127 
b; HD 216536 b) 


1. Introduction 


Searches for planets around stars that evolved off the Main Sequence (MS) represent 
an important ingredient of the exoplanet research. They extend the known population of 
planets to those orbiting massive and intermediate-mass stars. They also provide the best 
tool to study star-planet interactions induced by stellar evolution. Finally, these searches 
furnish high quality spectroscopic data that can be used in a wide range of studies related 
to stellar evolution. 


Intermediate-mass and massive MS stars have high effective temperatures and rotate 
rapidly. With only a few, rotation broadened spectral lines, these stars are not suitable 
for high precision radial velocity (RV) searches for planetary companions. Consequently, 
planets around intermediate-mass and massive MS stars discovered occasionally in transit 
searches (cf. |Borucki et al. 2011 Schwamb et al. 2013) are very difficult to confirm with 


RV measurements. The main source of data on such planetary systems are RV searches 
that focus on subgiant and giant stars which are evolving off the MS, cooling down, and 
considerably slowing their rotation. The most massive hosts of planetary systems come 
almost exclusively from such surveys (e.g. Sato et al. 2007, 2010, 2012, 2013). Surveys of 
this kind have been uncovering a population of planetary systems that are descendants of 
an undetermined original sample of planets around intermediate-mass and massive stars, 
carrying with them imprints of the effects of several Gyr of dynamical evolution enforced by 
an evolving star. They have demonstrated, for example, a paucity of planets within 0.5 AU 
of their parent stars ( Johnson et al. ||2007 Sato et aT7||2008 ), a borderline currently broken by 
Kepler 91b, a planet at only a/R* = 2.46 (Lillo-Box et al. 2014; Barclay et al.| 2014). 

The other frequently discussed effects are the higher frequency of planets around more 
massive stars and the stellar mass - planetary system mass relation (Lovis & Mayor 2007 


1 Based on observations obtained with the Hobby-Eberly Telescope, which is a joint project of the Uni¬ 
versity of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians- 
Universitt Mnchen, and Georg-August-Universitt Gttingen. 
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Johnson et al. 

2010, 

2011 

). Such stars also seem to prefer massive, m p > 1.5 Mj, planets 

(Sato et al. 2007) 

Niedzielski et al. 2009; 

Bowler et al. 

2010) and a more frequent occurrence 


of brown dwarf - mass companions, which are deficient around solar-mass stars on the MS 


Liu et al. 

2008 

Niedzielski et al. 

2009 

Mitchell et al. 

2013 

Niedzielski et al. 

2013) 


Disentangling information on the original planet population from the available data on 
planets around giants requires a detailed understanding of the interaction between an orbiting 
planet and an evolving star. In fact, importance of the influence of tidal interactions on 


orbital evolution has been recognized long before first exoplanets were discovered (see Kopal 


1972 Counselman 1973; Zahn 1977 and references therein). Livio (|1982) studied the net 


effect of accretion drag and mass loss in a star-planet system. The first attempts to predict 


planet survival upon stellar evolution were presented by Choi (1980); Choi & Vila (1981), 


Livio & Soker (1983, 1984) and Vila (1984). Detailed effects of brown dwarf or planet capture 


by an evolving star were numerically modeled in Siess & Livio (1999a b). The future of the 


Solar System in the context of the Sun’s evolution was studied by Duncan & Lissauer (1998); 


Rybicki & Denis 

(2001) 

Schroder & Connon Smith 

(2008 


A more general picture of stellar evolution induced star-planet interaction, including 


m 

Rasio et al. ( 

1996); 

Villaver & Livio 

(2007) 

Carlberg et 

Nordhaus et al. 

(2010 

). Three possible outcomes of orbital 


evolution of planets affected by tidal forces and stellar evolution, somewhat analogous to 
those in planet-satellite system ( Counselman|1973 ) were presented in Villaver & Livio (2009). 


Most recent studies include such effects as evolution of the primary’s spin (Nordhaus & 


Spiegel 2013) or interactions with the Galactic environment (Veras et al. 2014). The influence 
of stellar wind and tidal force prescriptions on planetary orbit evolution has been recently 


illustrated in Villaver et al. (2014). 


In this paper we present three planetary-mass companions around solar-mass stars 
from the ongoing Penn State-Toruii Planet Search (PTPS) with the Hobby-Eberly Tele¬ 


scope (Niedzielski et al. 2007; Niedzielski & Wolszczan 2008): BD+49 828, HD 95127 and 


HD 216536, and discuss their fate through their host evolution. We also summarize available 
data on planetary systems of evolved stars and briefly discuss properties of planetary systems 
at various stages of their hosts stellar evolution. 


The plan of the paper is as follows: in Section [2] we describe the observations and 
present our targets and their basic parameters, the measurements and the analysis of RVs. 
In Section [3] we present the results of the Keplerian modeling of RV variations, while other 
sources of stellar variability are considered in Section [4j Finally, in Section [5] we discuss 
our results, the fate of planets orbiting giant stars, and some additional problems related to 
evolved planetary systems. 
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2. Observations, stellar parameters and radial velocitiy measurements 

2.1. Observations 


Spectroscopic observations presented here were made with the Hobby-Eberly Telescope 
(HET) (Ramsey et al.j 1998) and its High Resolution Spectrograph (HRS) (Tull] 1998) in the 
queue scheduled mode (Shetrone et al. 2007). The spectrograph was used in the R=60,000 
resolution mode with a I 2 gas cell inserted into the optical path, and it was fed with a 2 arcsec 
fiber. The configuration and observing procedure employed in our program were practically 


identical to those described by Cochran et al. (2004). 


BD+49 828 (AG+49 353) is a K0 giant (Dieckvoss & Heckmann ||1997 ) with B = 9.43 ± 
0.03 (Hpg et al. 2000[ ) and V = 9.38 ± 0.02. We collected 42 epochs of precise RVs for this 
star over 3134 days (8.58 yr), between MJD 53225 and 56359. Typical signal-to-noise ratio 
(S/N) of about 200 was reached in 540-1600 s depending on atmospheric conditions. 


RVs of HD 95127 (BD+44 2038), a K0 giant ( 

Cannon & Pickering 

1918-1924 

) with 

B = 10.825 ± 0.045 and V = 8.15 ± 0.01 rfHpg et al. 

2000 

) were measured at 41 epochs over 


the period of 2929 days (8 yr) between MJD 53024 and 55953.The S/Ns of the combined 
star and iodine spectra ranged from 80 to 250. The S/N of the stellar template was 480. 
The exposure time ranged between 170 and 800 s, depending on the atmospheric conditions. 


Hipparcos parallax from van Leeuwen (2007) puts this star at a distance of 332 ± 110 pc. 


HD 216536 (BD+43 4329 ) is a K0 giant ( 

Cannon & Pickering 

1918-1924 

) with B = 

10.42 ± 0.03 and V = 9.23 ± 0.02 (Hpg et al. 

2000 

). The star was observed over 2756 days 


(7.54 yr), between MJD 53544 and 56300. The typical S/N of 250 was reached in 455 - 
1200 s exposures depending on actual conditions. 


2.2. Stellar parameters 


All the three giants discussed here belong to the PTPS Red Giant Clump sample. 


They have been studied in detail by Zielinski et al. (2012), who estimated their atmospheric 
parameters by means of the method of Takeda et al. (2005a|b) and intrinsic colors. In a 


more detailed chemical abundances analysis of Adarnow et al. (2014) with SME (Valenti & 


Piskunov 1996) no anomalies were detected and new rotation velocities presented. 


Luminosity of HD 95127 was calculated from the Hipparcos parallax and we found it 


to be surprisingly high for a log g = 3.82 ± 0.27 and T eS = 4807 ± 90 K star (Zielinski 


et al. 2012). To clarify that mismatch we reanalysed all available HET/HRS spectra with 
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no iodine imprinted and performed spectroscopic analysis with the method of fTakeda et ah 
(200 5a|b ) again. Resulting, significantly different atmospheric parameters are presented 

Table CD 


m 


Using these parameters the mass and age were derived with the Bayesian method based 


on 


Jprgensen & Lindegren (2005) formalism and modified by da Silva et al. (2006), in which 


we used theoretical isochrones of Bressan et al. (2012). For the other two stars luminosity 


was obtained, together with the mass and age, applying the Bayesian approach described 
above. 

Stellar radii were estimated either from the derived masses and logg values obtained 
from the spectroscopic analysis or using derived luminosities and effective temperatures. 


With the rotational velocities from Adamow et al. (2014) and stellar radii, masses and 


luminosities derived here we have obtained estimates of maximum rotation periods, as well 
as periods and amplitudes of p-mode oscillations from the scaling relations of Kjeldsen & 
Bedding (1995). All the information on targets is presented in Table [l} 


Positions of the three stars discussed in the paper on the Hertzsprung-Russel Diagram 
(HRD) are presented in Figure [I] together with the entire PTPS sample. In this Figure we 
show 332 stars from Zielinski et al. (2012), 368 stars from Niedzielski et al. (in prep.) and 
138 stars from Deka et al. (in prep.) for which spectroscopically determined atmospheric 
parameters as well as masses, radii and luminosities are available. 


2.3. Radial velocity measurements. 


The HRS is a general purpose spectrograph which is neither temperature nor pressure 
controlled and the calibration of RV measurements with this instrument is best accomplished 
with the I 2 cell technique. For the purpose of measuring the RVs and line bisector velocity 


span (BS), we have combined the method described by (Marcy & Butler 1992 Butler et al. 


1996) with that of cross-correlation used by (Qucloz 1995; Pepe et al. 2002). A detailed 


implementation of our approach is described in Nowak (2012) and Nowak et al. (2013). The 
typical precision of a few ms' 1 achieved that way made it quite sufficient to use the Stumpff 


(1980) algorithm to refer the measured RVs to the Solar System barycenter. 


The precision of both RVs and BSs as well as the long-term stability of our RV mea¬ 
surements has been verified by analyzing the data derived from monitoring of stars that do 
not exhibit detectable RV variations and stars for which the RV or BS variations have been 


well established. One example of such a calibrator, discussed in Niedzielski et al. (2009), 
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is the KO giant BD+70 1068 , which has (Try = 12ms -1 consisting of the intrinsic RV un¬ 
certainty of 7ms -1 and the approximate 10ms -1 RV amplitude variation due to solar-type 
oscillations Kjeldsen & Bedding (1995). Nowak et ah (2013) have carried out a more detailed 
analysis of our RV precision and stability using the RV measurements and orbital solutions 
for HD 209458, HD 88133. The precision of our BS measurements has been demonstrated 
for HD 166435 by the same authors. 


3. Keplerian analysis 


Orbital parameters for the companions to the three stars were derived using a hybrid 
approach (e.g. (Gozdziewski et al. 2003; Gozdziewski & Migaszewski 2006[ [Gozdziewski et 


al. 2007)), which combines the PIKAIA-based, global genetic algorithm (GA; Charbonneau 


(1995)) with a faster and more precise local method. Given a periodic signal found in the RV 
data with the Lomb-Scargle (LS) periodogram ( Lomb||1976 Scargle |l982 ; Press et al. 1992), 
we begin with launching a search for possible orbital solutions over a wide range of parameters 
using the GA code. The GA semi-global search usually establishes a narrow parameter range 


that defines the search space for the next step, in which the MPFit algorithm (Markwardt 


2009) is used to find the best-fit Keplerian orbit delivered by RVLIN (Wright & Howard 


2009) modified to allow the stellar jitter to be fitted as a free parameter (Ford & Gregory 


2007 Johnson et al. 2011). The RV bootstrapping method (Murdoch et al. 1993 Kuerster 


et al. 1997 Marcy et al. 2005; Wright et al. 2007) is employed to assess the uncertainties of 
the best-fit orbital parameters defined as the width of the resulting distribution of 10 6 trials 
of scrambled data between the 15.87th and 84.13th percentile. 

The false alarm probability (FAP) of the final orbital solution was derived by repeating 
the whole hybrid analysis on 10 5 sets of scrambled data. 


3.1. BD+49 828 

The RV measurements for this star are presented in Table [2] and Figure [2j The RV 
variability is characterized by an amplitude of 73 ms -1 and mean uncertainty of 8.5 ms -1 . 
The LS periodogram of these data (Figure [8]) reveals only one significant peak at a period 
of ~ 2410 days at the false alarm probability FAP ^5 < 0.01. 

The final fit of a Keplerian orbit to the data resulted in the post-fit rms residua of 
11.6 ms -1 and jitter of 4.4 m s -1 in very good agreement with our estimates. The fit produced 
the solution with y 2 = 1.35, and, after 10 5 trials the FAP < 10 -5 . The final best-fit Keplerian 
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parameters are listed in Table [5] and the model orbit is plotted together with the data in 
Figure [5j No additional periodic signal was detected in the post-fit RV residua (Figure [ 8 ]). 

This solution indicates the presence of a rn-i sin i = I.GIq^Mj body in a 4 . 2 +q; 2 2 AU, 
moderately eccentric, e = 0.35 ^q; 2 q orbit with the period of P = 259CRfgo days. Obviously, 
the long orbital period of BD+49 828 b is currently not very well constrained. Clearly, no 
less than 8 years of additional observations are needed to determine it more precisely. 


3.2. HD 95127 

RVs of HD 95127 (Tableland Figure[3]), with an amplitude of 404.8 m s -1 and estimated 
mean RV uncertainty of 5.98 ms -1 , reveal a clear 483-day periodicity (FAP is <0.001) in the 
corresponding LS periodogram (Figure [9]). 

The orbital solution includes a planet with the minimum mass of m 2 sini= 5 . 0 l!o 44 Mj, 
in a 482^5 day, 1.28lom AU, e = 0.1llo!o6 or bit around the star with the y 2 = 1.14 and 
FAP< 10 -5 (Table [HJ Figure [ 6 ]). The post-fit residua show no significant periodicities in the 
corresponding LS periodogram (Figure [9]). 

We also note that the 47.5 m s -1 jitter and 50.9 m s -1 post-fit residua from the final model 
are consistent with the stellar jitter interpreted as under-sampled p-mode oscillations (highly 
skewed distribution of jitter in our bootstraping analysis prevented a reliable uncertainty 
estimate). 


3.3. HD 216536 


The RVs for HD 216536 are presented in Table [4] and plotted in Figure |4j In this 
case, the mean uncertainty and amplitude characterizing the RV variability are 6.4 ms -1 
and 144 ms -1 , respectively. The LS periodogram of the data (Figure 10) shows a single 
significant period of 148 days with the FAP ^5 <0.001. 


The hybrid modelling of these data, resulted in Keplerian fit with the post-fit rms 
residua of 23 ms -1 and a jitter of 17.9 ms -1 (again presenting a highly skewed distribution), 
consistent with our estimates of the amplitude of unresolved p-mode oscillations based on 
stellar mass and luminosity. The final, best-fit model (y 2 = 1.29, FAP< 10 -5 ) parameters 
are listed in Table [5] and the model RV curve is plotted in Figure [7} No other periodic signal 
is present in the RV residua. 

A minimum mass of the companion to HD 216536 is m 2 sin i = 1 . 47 ^q' 2 2 Mj, and its 



close-in, 0.609 ^o;qo2 AU, P = 148. 6 / 0)7 day orbit has a moderate eccentricity, e = 0.38 /o/q. 


4. Stellar activity analysis: line bisectors, Ha variations and photometry. 


Both RV searches for planets and studies of star-planet interaction are affected by the 
intrinsic variability of host stars causing RV or BS variations not related to stellar reflex 
motion caused by the presence of an orbiting planet. Changes in spectral line shapes arising 
from motions in the stellar atmosphere, related to nonradial pulsations or inhomogeneous 
convection and/or spots combined with rotation can mimic low-level long period RV vari¬ 


ations. Significant variability of red giants has been noted already by Payne-Gaposchkin 


(1954) and Walker et al. (1989) and made the nature of these variations a topic of numerous 


studies. 


Hatzes & Cochran (1993) showed that the low-amplitude, long-period RV variations 


may be attributed to pulsations, stellar activity - a spot rotating with a star, or low-mass 
companions while the presence of short-period, possibly p-mode oscillations as demonstrated 


m 


Hatzes & Cochran (1994). 


The solar-type, p-mode oscillations are easily observable in high precision, photometric 


time-series measurements and they have been intensely studied based on the COROT (Baglin 


et al. 

2006 

) and KEPLER ( 

Gilliland et al. 

2010 

) data, leading up to precise stellar mass 


determinations (De Ridder et al. 2009 Bedding et al. 2010; Kallinger et al. 2010 Hekker et 


al. 2011). With the typical timescales for giants, ranging from hours to days, such oscillations 


typically remain unresolved in low-cadence observations focused on long-term RV variations 
and they contribute as an additional uncertainty to the RV measurements. 

Long-period RV variations in some giants have been successfully demonstrated to be 


due to the presence of low-mass companions ( 

Campbell et al.|1988| ' 

hatzes et al. 2003 

Hatzes 

& Cochran 

1993 

Hatzes et al. 

2006 

Vogt et al. 

2000 

Butler et al. 

2001 

Frink et al. 

2002 ). 


According to the exoplanets.eu website, 55 giant stars with substcllar-mass companions have 
been identified so far. However, in some giants, the nature of the observed RV variations 


indicating the presence of long-term, secondary periods (Payne-Gaposchkin 1954 Houk 1963) 


remains a riddle (Nicholls et al. 2009). 


Therefore, a reliable planet detection requires all the known sources of the observed RV 
variations to be investigated in a considerable detail. We note, however, that in the case 
of the three stars discussed in this paper, substantial eccentricities of the best-fit Keplerian 
orbits alone make stellar pulsations based alternative explanations less likely. 
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Ca II H&K resonance lines, if measurable, are good indicators of stellar activity. Un¬ 
fortunately, in the HET/HRS spectra of red giants these lines are too weak to be usable. 
Therefore, we have developed a cross-correlation technique to monitor stellar activity with 
BS measurements derived from exactly the same spectra that have been used for RV deter¬ 


minations (Nowak 2012 Nowak et al. 2013). 


As an independent stellar activity indicator we also measured the Ha index (/h q ) from 


HET/HRS red spectra based on the approach presented by Gomes da Silva et al. (2012) and 


Robertson et al. 

(2013 

and references therein), following the procedure described in detail 

in 

Maciejewski et al. 

(2013 

)• 


Both activity indicators, BS and Jh q were computed at the same epochs and from the 
same spectra as the RVs but such subtle effects as instrumental profile variations were not 
accounted for in these measurements. Consequently, they are, to some extend, subject to 
instrumental effects presumably not present in the RV measurements. As the Ha line is also 
rather weak in K giants it is especially prone to instrumental effects. To monitor such effects 
we also measured the index of the Fe I 654.6239 nm line (Jpe), which is insensitive to stellar 
activity. 


Also, the PTPS stars are relatively bright, which makes their photometric measure¬ 
ments available from the published surveys. These three kinds of data provide most of the 
information that we use to diagnose any possible activity of our target stars that may mimic 
RV variations due to the Keplerian motion. 


4.1. BD+49 828 


The maximum rotation period of BD+49 828, given the estimated projected rotation 
velocity of u rot sin R = 1.7 ± 0.8kms _1 and stellar radius of 7.6 ± 1.3 R 0 , is 226 ± 114 days, 
much less than the orbital period, which makes rotation of a spot on the stellar surface an 
unlikely scenario. 


There are two archival photometric time-series available for this star. The Northern Sky 
Variablity Survey (NSVS; Wozniak et al. 2004) contains a set of 81 epochs of photometric 
measurements made over 216 days between MJD 51414 and 51630, which give the mean 
brightness and rms scatter of 9.23 mag and 0.028 mag, respectively. The Wide Angle Search 
for Planets (WASP; Butters et al. 2010) archive contains a much larger data set of 5557 
epochs of photometric observations made over 1262 days between MJD 53196 and 54458. 
These measurements are partially contemporaneous with our RV measurements and are 
characterized by the respective mean brightness and the rms scatter of 9.659 mag and 0.019 
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mag. None of these data series contain a periodic signal close to the best-fit Keplerian orbital 
period. The LS periodogram of the WASP observations is shown in Figure [8} An absence 
of any periodic signal in the existing photometric time series suggests that pulsations or a 
spot rotating with the star cannot be the source of the observed RV variations. However, 
the time covered by photometric observations is short compared to orbital period. Therefore 
we conservatively assume that all the observed photometric scatter originates from a spot 
of a size of / = 1.9% of the rotating stellar surface. Then, according to Hatzes (2002), we 
might expect a periodic variation of the RVs at a full amplitude of 64ms -1 . This would be 
similar to the observed RV variations and it would generate BS variations at the amplitude 
of 45 ms -1 correlated with RV. To resolve such a situation one needs to carry out a detailed 
study of BSs and search their variations for periodicities that would correlate with those 
observed in the RV measurements. 


The BS measurements for this star have a mean value of 18±33ms _1 and amplitude of 
142 m s" 1 with the mean uncertainty of 19 m s -1 . The Parsons correlation coefficient between 
the RV and the BS variations amounts to -0.09 indicating that they are independent. In 
addition, the BS variations do not exhibit any measurable periodicities (the LS periodogram 
of the BS variations is shown in Figure [8]) Finally, the 2590+igg -day orbital period of 
BD+49 282 b is much too long to be related to the observed rotational period of the star, 
which is more than 11 times shorter than that. 


Our Jhq, measurements for BD+49 828 present a mean value of 0.0365 and range between 
0.0381 and 0.0353, i.e. explicit variations of 7.5%. At the same time Ip e has a mean value 
of 0.01147 and ranges from 0.01169 to 0.01132 i.e. varies by 3.3% accounting for nearly a 
half of the observed Iu a variations. The LS periodograms of both indices are shown in Fig. 
[8j The In a shows a significant periodicity of about 1750 days (~ 5 yr) and Ip e data show 
a signficant period of about 1 yr, apparently of seasonal origin. The long period signal is 
well separated from the RV period of 2590^^gg days. The 1 yr signal is, however, close to 
the estimated rotation period of the star. Both LS of line profile indices show a similar 
pattern suggesting strong influence of instrumental effects enhanced in the case of J Hq , by 
stronger signal-to-noise ratio. Indeed, the Parson’s correlation coefficient between RV and 
/hc is r=-0.344 only while the correlation between Jh q and Jf c is r=0.620, much above the 
critical value of 0.408 (confidence level of 0.01). We may safely assume therefore that the 
observed /h„ variability is induced by instrumental effects and carries no evidence of stellar 
activity. 

Based on the discussion above, we conclude that the RV variations in BD+49 828 cannot 
be explained in terms of rotation induced stellar activity or pulsations. Consequently, a 
Keplerian modeling of the measured RV variations of the star is justified and the presence of 
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an orbiting low-mass companion represents the most natural interpretation of the data. The 
observed BS variations must originate from another source which is either random or not 
resolved with our observations. The origin of Iu a and Jp e variations seems to be instrumental 
but it’s nature is not clear. 


4.2. HD 95127 


With the estimated projected rotation velocity of i> rot sini* = 2.6±0.7kms _1 and radius 
of 20 ± 9 R 0 the maximum rotation period of HD 95127 is 389 ± 204 days, much less that 
the orbital period. 


If the observed RV variations were indeed a result of a spot rotating with the star the 


observer RV semi-amplitude would require a spot of 9% of the stellar surface (Hatzes 2002). 


Such a spot would be easily detectable with moderate quality photometry. The only existing 


extensive V-band photometric monitoring for HD 95127 comes from SWASP (Pollacco et 


al. 2006) where 540 epochs are available over only about 100 days between MJD 54091.1 


and 54194.9. The mean brightness of the star was 8.394 ± 0.012 mag and the data show an 
amplitude of 0.033 mag only, three times less than expected. 95 Hipparcos (ESA 1997) H p 
measurements taken over three years between MJD 47877 and 48961 show also a constant 
brightness of H p = 8.307 ± 0.014 mag and no significant periodicity (Fig. [9]). We can 
conclude therefore that the existing photometric data do not support the spot hypothesis. 
The lack of any significant periodic signals in photometric data excludes also pulsations as 
a possible source of RV variations. 


The spot hypothesis predicts BS variations of semi-amplitude of 23.7ms -1 (Hatzes 
2002) while the observed BS semi-amplitude is 42ms -1 . However, the uncertainties in BS 
are much larger than in RV (16.6 m s -1 on average) and the mean value of BS is 23±23 m s _1 . 
The LS analysis of BS data shows no trace of periodicity at the rotation period or at the 
Keplerian best fit period (Fig. [9]). Moreover, the Parson’s correlation coefficient between 
RV and BS is r=0.25 with the critical value of 0.4, and no relation between these quantities 
can be stated. Again, the BS data do not support the spot hypothesis. 


The Jh q index for HD 95127 has a mean value of 0.0339, and varies by 8.5%. At the 
same time Jp e index presents a mean value of 0.00983 and varies only by 1.9%. None of 
indices shows any significant periodic signal (Fig. [9]). Also the correlation analysis carries 
no useful information on activity as both correlations between RV and Jh q and between Jh» 
and Jpe, r=0.419 and r=0.424, respectively are at the level of critical value (0.403). 


Therefore we conclude that the only significant periodic signal present in the RV data 
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is due to an orbiting planet and there is no indication of pulsations or spot rotating with the 
star that can be inferred from our data. 


4.3. HD 216536 


The estimated projected rotation velocity of v rol sin % = 2.6±0.5kms~ 1 and the radius 
of 12.5 ± 4.0 R 0 of this star result in the maximum rotation period of 243 ± 91 days. HD 
216536 had been observed by the NSVS (Wozniak et al. 2004) 5.6 years before the beginning 
of onr measurements. Photometric observations of the star at 90 epochs between MJD 
51304 and 51483 have been collected. These measurements are characterized by the mean 
brightness of 9.118 mag and the rms scatter of 0.029 mag. As shown in the LS periodogram 
of Figure [TOj no significant periodicity is present in the data. The observed photometric 
scatter could be explained in terms of a f= 2.9% spot rotating with the star. This would 
result in periodic, 185ms -1 and 25ms -1 variations in the RVs and BSs, respectively (Hatzes 


2002). The fact that the predicted, spot generated RV variations are similar in amplitude 


to those actually observed requires a careful analysis of the BS behavior. 


The measured BSs exhibit variations that have the mean value of 14 ± 22ms -1 , the 
amplitude of 106ms _1 , and the mean uncertainty of 14ms' 1 . The Parsons correlation 
coefficient of the RV and BS variations amounts to 0.29 which means that they are not likely 
to be correlated. 


In addition, BS variations do not show any measurable periodicities. Consequently, the 
absence of any correlation between RVs and BSs and lack of periodic variations in both BS 
and existing photometry invalidates the spot hypothesis. 

Our measurements of I\] n index for HD 216536 reveal a mean value of 0.0354, and 
amplitude of variations of only 0.0021 (6%). The Ip e variations are even weaker - 2.6%. 
Both indices show very similar pattern of LS and there is a significant signal at 984 days 
in Jh q data. The signal is well separated from the Keplerian period and the estimated 
maximum rotation period of 347±115 days. The Parson’s correlation between /n a and Jpe, 
r=0.595 is above critical (r=0.376) which proves the instrumental nature of the signal. The 
RV variations are unrelated to that effect as the correlation coefficient between RV and Iu Q 
is only r=0.055. 

We conclude that, according to the available evidence, the RV variations in HD 216536 
cannot be related to rotation induced stellar activity or stellar pulsations. Consequently, it 
is the most logical to assume that these variations are due to the orbital motion of a planet- 
mass companion around the star. As for the other two stars, the observed BS variations 
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must have another, unrelated source which is either random or not resolved in time with our 
observations. 


5. Discussion. Planetary systems of evolved stars. 


The three new stars with planets presented here, provide a representative sample of the 
wealth of planetary system architectures observed in the evolved systems and illustrate how 
searches for planets around such stars can address various aspects of stellar and planetary 
systems evolution. 


By exploiting the approach of Ford & Gregory (2007) and Johnson et al. (2011) in 


Keplerian orbit fitting we obtained independent estimates of stellar jitter, whose amount is 
consistent with the amplitudes of p-nrode oscillations calculated from the scaling relations of 


Kjeldsen & Bedding (1995). This proves that the L/M ratios for the target stars and hence 
their positions on the HRD are correct. 


HD 95127, a bright, logL/L 0 = 2.28 ± 0.38, solar-mass M = 1.2 ± 0.22 M 0 giant, 
with R = 20 ± 9R 0 and [Fe/H]= —0.18 ± 0.05 hosts a mshii = S.OI^q^Mj planet in 
a 1.28 ^qq4 AU, e = O.ll^QQg orbit. The relatively low mass of the star and its expected 
large radius at the tip of the Red Giant Branch (RGB) set the current orbit of its planetary 


companion well within the engulfment zone (Villaver & Livio 2009 Kunitomo et al. 2011 


Villaver et al. 2014). This presents an independent proof that the star is still evolving up 


the RGB. 

HD 216536, a M/ M 0 = 1.36±0.38 and [Fe/H]= —0.17±0.09 star, within the parameter 
uncertainties, is also evolving up the RGB. Most likely this star is still before the helium 
ignition. It harbors one of the most close-in planets orbiting a giant, which supports our 
conclusion concerning its evolutionary stage. Given its 0.6091 q;qq 2 AU orbit, this 1.471 q 0 2 Mj 
planet is one of the very few warm Jupiters around giants, prone to strong tidal interactions 


and most likely to be ingested before its host reaches the tip of the RGB (Villaver & Livio 


2009 

Kunitomo et al. 

2011 

Villaver et al. 

2014) 


The position of BD+49 828 on the HRD suggests that this M/ M 0 = 1.52±0.22, slightly 
metal deficient ([Fe/H]=-0.19±0.06) star, is evolving up the RGB, possibly undergoing the 
first dredge-up. Its m sin i = I.O^q^Mj companion, 4.2 +q; 2 2 AU away of the star, is one of 
the most distant planets orbiting giants detected so far. It is very likely that the planet 


will not be affected by stellar evolution and similar to HD 4732 c (Sato et al. 2013) or HD 


120084 b (Sato et al. 2013). Taking into account the stars’s uncertain luminosity estimate, 


the planet may also be located within the habitable zone as defined in Kopparapu et ah 
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(2013) 


To put the three planetary candidates discussed in a wider perspective, we plotted them 


as red circles in the log(g) versus log(a/.R*) plane in Figure 11 Because the orbital radius, 
a/R*, is a major factor scaling the strength of tidal star-planet interactions, this plot may be 
regarded as an illustration of the tidal evolution of planetary systems. In the same Figure we 
have also placed all known exoplanet^] and KOIs that orbit stars with masses in the range 
of 1-2 M 0 as black circles. Finally the total sample of planetary companions, for which with 
at least semi-major axes, and stellar masses and radii are available is shown as grey circles. 
For all stars, including the three ones presented here, uniform logo’s were calculated from 
stellar masses and radii adopted in the discovery papers. We note, that for some planetary 
systems hosts the spectroscopic logg do not match those calculated from stellar mass and 


radius. In particular this concerns 7 CMa (Wittenmyer et al. 2011), another long-period 
planet for which the host star log gM,R = 3.9. 

Also included in this Figure is the minimum distance to avoid engulfment for a 1 Mj 


planet orbiting a 1.5 M 0 star under the Schroder & Cuntz (2005) mass-loss prescription (red 
line) and the minimal orbit beyond which the planet is not affected by the tidal forces (blue 


line) from Villaver et al. (2014). The dashed line shows the a/i?* = 3 border inside of which 


tidal interactions are expected to lead quickly to ingestion. Note that the minimum distance 


to avoid engulfment for a lMj planet orbiting a 1.5M 0 star from Villaver et al. (2014) 


(red line) fits the general slope of decreasing a/R with log g for known planetary systems in 
Figure [Tl] In fact, up to a very late stellar evolution phase at RGB where log g < 1 and tidal 
ingestion is happening this relation simply reflects the steady growth of a stellar radius. 


An interesting feature of Figure [Tl] is that after the discoveries of HD 102956 b (Johnson 
et al. 2010), Kepler-56 b, c ( jSteffen et al.| 2013), Kepler-391 b, c (Rowe et al. 2014), KOI- 
1299 b ( Ciceri et al. ||2014 ; Ortiz et al. 2014) and HIP 67851 b ( |Jones et al.||2014 ), the gap in 
orbital separations of planetary systems orbiting subgiants (here: log <7 = 3.5 ± 0.5) pointed 


out in Nowak et al. (2013) no longer exists. A relatively low number of systems around that 
kind of stars is obviously a consequence of a low population of that part of the HRD (see for 


example Figure 3.5.5 in ESA 1997) 


The ingestion itself is supposed to be a quick process as clearly illustrated in Figure 


11 Only three dwarfs (here: log g = 4.5 ± 0.5) (Wasp-12 b (Hebb et al. 2009), Kepler-78 b 
( Sanchis-Ojeda et al.||2013 ), Wasp-103 b ( Gillon et al.||2014 )) out of 1037 ones or ~ 0.3% are 
located in that range. The same is true for only one (Kepler-91 b) out of 44 or ps 2% giants 
(here: log g = 2.5 ± 0.5). 


2 All data from exoplanets.eu 
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The most evolved planetary systems with logg < 1.5, very close to engulfment, are: 
HD 220074 b, HD 208527 b ( [Lee et al.|[M3| ), [3 Cnc b ( [Lee et al.|[2014] ), HD 96127 b 


(Gettel et al. 2012), (3 Umi b (Lee et al. 2014), all within a/R * = 9. Interestingly enough, 
except for Kepler-91 b there is still no planet orbiting an evolved star with logg < 3 within 
a/R+ 7 (the tightest orbit being HD 220074 b at a/R* = 6.93). This may suggest 
that tidal interactions between most extended evolved stars and their planets are more far- 
reaching. 


It is important to note, that the part of Figure 11 occupied by bright giants, (here: 


logg = 1.5 ± 0.5) may be populated by both planets hosted by RG stars before RGB tip 
and engulfment, and post RGB-tip stars planets. Although the latter planets have avoided 


tidal ingestion at the RGB tip, their orbits must have decayed (channel (2) in Villaver 
& Livio 2009, see discussion of their Figure 1), in some cases to well within 1 AU, and 


may subsequently be subjected to ingestion during their hosts Asymptotic Giant Branch 
evolution. Undoubtedly, all these stars deserve more detailed spectroscopic analysis to asses 
their evolutionary stages and hence the dynamical history of their planets in more detail. 

All planets below the red line with at least one Jupiter mass are expected to be eventually 
ingested by their hosts before they reach the tip of the RGB. An obvious conotation from 
Figure [Tl] is that vast majority of known planetary systems is located within the minimum 
distance to avoid engulfment (red line) and will not survive RGB evolution of their hosts, 


in agreement with Carlberg et al. (2009) and Nordhaus et al. (2010). Only very few known 


planetary systems contain planets in orbits, in which are not affected by the tidal forces 
(above the blue line). This is 105 out of 1205 planets or ~ 9% of planets in the exoplanet.eu 
catalog, for which at least semi-major axes, M* and R* are available. Most of the safe planets 
orbit dwarfs (77 out of 1037 or ~ 7%). A few such systems exist around subgiants (10 out of 
99 or « 10%) with the most evolutionary advanced of them being HR 228 b and c (Sato et al. 


2013). Only 5 (out of 44 or « 11%) of such systems are present around giants: HD 139357 b 


( Dollinger et al.||2009 ); u Oph c and omi UMa b ( Sato et ah|2012 ); HD 120084 b (Sato et al. 


2013); HD 14067 b (Wang et al. 2014). BD+49 828 b presented here belongs to that group 


as well. 

The obvious shortage of known safe planetary systems is most likely a consequence of 
current observational limitations. It illustrates the importance of long-term planet search 


projects like, for example, PTPS, the Okayama Planet Search (Sato et al. 2005) or new 


projects like Friends of Hot Jupiters (Knutson et al. 2014). It also has serious consequences 


for identification of asymmetric planetary nebular formation processes (Harpaz & Soker 1994 


Soker 

1996 

Livio & Soker 

2002) 


The giant planets are formed outside the snow-line, at a distance of 6 - 12 AU in the 
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m 


case of 1 — 2 M 0 stars ( Kennedy &; Kenyon||2008 ). They migrate inward (Armitagc 2007) 
the protoplanetary disks through type II migration but other migration scenarios: planet- 


planet scattering (' 

Jasio & Ford 

1996 

Wu & Murray 

2001 

5) are also possible. 


Ford et ah 2001) or Kozai migration (Kozai 1962 


the protoplanetary disk is dispersed. The efficiency of the migration process, and especially 
a relative efficiency of various mechanisms is basically unknown due to the lack of known 
long-period planets. We note the importance of Rositer-McLaughlin effect studies in that 
context (Simpson et ah 2010). 


An inventory of planetary system at various evolutionary stages of their hosts is pre¬ 
sented in Table [6j in which four sets of data are shown: all stars for which all parameters 
listed in Table [6] are available, and a subset of these stars limited to 1 — 2 M 0 , as well as sub¬ 
sets of those two sets with the RV detected planets only. In case of dwarfs, we see the impact 
of short-period transit planets as a huge variation between the mean value and the median 
of minimum planetary semi-major axis, hence the need for a ,,RV limited” subsamples. 


Several interesting features of the known sample of exoplanets emerge. In RV limited 
subsamples not much variation in semi-major amplitude or eccentricity is seen between 
planetary systems at various evolutionary stages of their hosts. Both a and e vary by 
about la between dwarfs and bright giants ruling out tidal circularization during the stellar 
evolution. 


Another important general feature of the known planet sample is the steady decrease 
of metallicity as we move from dwarfs to several billion years older giants. That reflects the 


changes in stellar populations we are dealing with, varying star formation rate (Tremblay 


et al. 2014) possibly reflected in the planet formation rate (Boss 2002 Johnson & Li 2012 


2013 Hasegawa & Hirashita 2014). 


As we move from dwarfs to more evolved hosts we deal with more massive stars on 
average as well. An average dwarf with a planetary system is a ~ 1M 0 (F or G spectral 
type) star while an average subgiant is already a ~ 1.5 M 0 (MS spectral type A-F) star, 
and a giant (MS A-type star) may be almost twice as massive as a dwarf (5-10cr difference). 
Stellar mass is not expected to increase during MS and RSG evolution and it is clear that we 
must be facing a strong selection effect here. Even in a narrow 1 — 2 M 0 mass range we still 
see more massive, evolved stars. On one hand, masses of single evolved stars are not easy 


to estimate, they may be simply overestimated (Mortier et al. 2013), which would be the 


easiest explanation of the suspicious lack of solar-mass evolved planetary systems hosts. On 
the other hand, taking into account that a solar mass star is expected to lose ~ 0.2 — 0.4 M 0 


in the RG stage ( 

Christy 

1966] 

Fusi Pecci et al. 

1993; 

Lebzelter & Wood 

2005 

Kafirai et al. 

2007 

McDonald et al. 

2011 

) and remembering that the fraction of RGB mass loss increases 
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with stellar mass (Kalirai et al. 2008), we see that the known evolved planetary system hosts 
might have been originally fa 2 — 3 M 0 MS stars, at the low - intermediate mass border or 
above. This points to their very uncertain, possibly quite complicated evolutionary history. 

The planetary mass increase for more evolved hosts seems to be the most prominent 
feature present in Table [6] An average giant hosts a companion about twice as massive 
than in the case of a dwarf, and a bright giant’s companion is 3 times more massive, on 
average. Certainly, that is not an observational bias caused by RV jitter present in giants 
as companions below Jupiter mass can be easily detected around those stars (for instance 


BD+48 738 b, Gettel et ah 2012). Assuming an average value of sin companions to bright 


giants are, on average, at the brown dwarf - planet borderline (see also Mitchell et ah 2013). 
ft is rather obvious that such an amount of mass cannot be accreted by a giant planet during 


its host’s RGB evolution (Duncan & Lissauer 1998). 


We note, however a lack of correlation between current host mass and its observed 
planetary system mass (Table [7]). ft does not necessarily mean that there is no relation 


between hosts and their planetary system masses postulated in Lovis & Mayor (2007). We 


may speculate that, to show such a relation, one should use initial masses of the planetary 
system hosts, instead of the current ones, corrected for mass-loss over their entire evolution, 
and account for the planetary system mass - loss through planet ingestion or ejection. If we 
assume, for simplicity, that an average subgiant from Table [6] (RV-limited, complete sample 
mean values of mass) lost 5% of its initial mass, whereas a giant and a bright giant would lose 
15% and 35%, respectively, we obtain a perfect host mass - planetary mass relation with a 
correlation coefficient r=0.95, which is statistically significant at 99% confidence level. This 
result is promising enough to justify a more detailed study. 
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Table 1. Stellar parameters of program stars. 


Parameter 

BD+49 828 

HD 95127 

HD 216536 

V [mag] 

9.38 ±0.02 

8.15 ±0.01 

9.23 ±0.02 

Spectral type 

K0 

K0 

K0 

7r[mas] 

— 

3.06 ±0.99 

— 

Teff [K] 

4943 ± 30 

4218 ± 69 

4639 ± 45 

log g 

2.85 ±0.09 

1.78 ±0.3 

2.36 ±0.21 

[Fe/H] 

-0.19 ±0.06 

-0.18 ±0.05 

-0.17 ±0.09 

v ro t sin i* [kms -1 ] 

1.7 ±0.8 

2.6 ±0.7 

2.6 ±0.5 

log L+/L q 

1.47 ±0.13 

2.28 ±0.38 

1.80 ±0.21 

M*/M 0 

1.52 ±0.22 

1.20 ±0.22 

1.36 ±0.38 

R-k/ R© 

7.6 ±1.3 

20 ±9 

12.5 ±4.0 

log (Age) [yr] 

9.37 ± 18 

9.74 ±0.27 

9.58 ±0.33 

P(sini*) _1 [days] 

226 ± 114 

389 ± 204 

243 ± 91 

Rose [ m s -1 ] 

4 c;+ 3 - 1 

37-25 

n+is 

ii -6 

Pose [days] 

n i q+o. 08 
u - ±o -0.05 


n qq+o.56 

U.Oa_o 25 



Tabic 2. Relative radial velocities and bisector span of BD+49 828 


Epoch (MJD) 

RV (m s' 1 ) 

0 -rv (m s x ) 

BS (m s x ) 

obs (m s *) 

53225.44499 

-20.0 

44.1 

79.1 

10.8 

53299.22458 

-2.9 

15.2 

-9.8 

26.6 

53593.40169 

4.0 

14.1 

-39.1 

24.1 

53608.39577 

0.5 

13.5 

-42.3 

21.6 

53610.36711 

15.2 

13.6 

-38.8 

19.5 

54071.31499 

17.8 

13.5 

-12.2 

15.9 

54087.06326 

18.9 

13.5 

-1.4 

16.6 

54101.25992 

17.0 

13.7 

34.6 

16.1 

54114.21772 

-5.8 

13.4 

10.6 

18.4 

54130.16924 

24.2 

13.6 

90.2 

15.0 

54148.11459 

-3.4 

13.3 

40.7 

17.4 

54449.07380 

4.6 

13.7 

7.6 

22.5 

54457.27139 

22.4 

13.8 

13.3 

20.3 

54519.10534 

34.1 

13.7 

58.9 

17.5 

54731.27897 

-3.5 

14.4 

22.9 

23.7 

54756.21112 

10.5 

14.2 

12.4 

22.3 

54781.17038 

19.5 

14.6 

34.9 

19.1 

55048.43124 

-27.9 

14.3 

23.6 

10.3 

55076.35465 

0.2 

15.3 

-51.9 

10.6 

55079.34718 

-9.3 

13.5 

-6.3 

10.6 

55107.26189 

8.7 

14.0 

15.0 

10.7 

55110.26035 

4.0 

13.7 

6.3 

10.9 

55179.07981 

15.4 

14.5 

17.9 

22.3 

55206.22428 

-25.1 

13.7 

4.2 

15.1 

55245.12414 

2.8 

14.6 

74.9 

22.4 

55247.11719 

-9.5 

14.1 

70.3 

21.0 

55443.34090 

-35.4 

15.0 

7.6 

32.0 

55478.25636 

-39.1 

13.7 

41.8 

16.5 

55501.18840 

-26.4 

13.8 

39.8 

19.6 

55518.16037 

-5.7 

13.5 

0.6 

17.0 

55775.45399 

-11.8 

13.1 

16.5 

16.9 



Tabic 2—Continued 


Epoch (MJD) 

RV (m s" 1 ) 

(Try (m s x ) 

BS (m s x ) 

o-rs (m s *) 

55791.37892 

-10.6 

13.5 

37.5 

19.9 

55799.38224 

-18.2 

13.1 

10.7 

17.8 

55825.30463 

-9.6 

13.6 

-3.5 

15.6 

55859.21333 

2.3 

13.7 

14.8 

13.2 

55931.23168 

-10.6 

13.6 

31.3 

13.7 

55973.12528 

5.1 

13.3 

40.2 

15.8 

56214.24030 

7.0 

13.6 

7.7 

20.6 

56257.13616 

3.9 

13.6 

17.7 

15.1 

56324.16275 

31.9 

18.5 

-35.1 

51.9 

56347.08323 

4.5 

13.6 

54.5 

18.0 

56359.08653 

7.8 

14.2 

46.3 

23.6 



Table 3. Relative radial velocities and bisector span of HD 95127. 


Epoch (MJD) 

RV (m s” 1 ) 

cr RV (m s x ) 

BS (m s x ) 

OBS (m s x ) 

53024.52304 

38.3 

7.1 

23.1 

17.8 

53370.35059 

-82.1 

9.1 

12.3 

24.0 

53867.22760 

-138.3 

5.5 

34.6 

13.3 

54048.49149 

47.8 

5.7 

5.4 

18.4 

54129.28005 

137.1 

6.6 

19.8 

17.9 

54153.43199 

120.2 

6.1 

-4.5 

16.9 

54165.40075 

32.2 

5.7 

58.9 

15.8 

54198.31076 

37.0 

4.3 

32.6 

13.3 

54212.26155 

91.8 

5.9 

16.1 

19.7 

54214.26450 

87.1 

4.8 

2.9 

18.3 

54233.22538 

-35.9 

4.6 

21.3 

14.0 

54249.17289 

-43.3 

5.2 

6.8 

15.9 

54249.17781 

-46.4 

4.3 

-8.3 

12.1 

54413.47777 

-88.4 

5.3 

55.3 

17.1 

54439.43738 

-51.3 

6.4 

57.8 

18.7 

54498.47970 

11.5 

7.1 

8.8 

21.1 

54532.18166 

107.8 

7.1 

55.3 

18.4 

54567.33124 

150.4 

6.9 

22.1 

18.5 

54568.09186 

99.1 

6.8 

63.9 

18.5 

54580.27156 

75.8 

5.8 

31.9 

15.6 

54618.14668 

112.7 

6.4 

41.6 

17.5 

54811.41721 

-40.9 

6.4 

-19.2 

19.5 

54842.33315 

-35.1 

5.9 

17.0 

18.4 

54867.24900 

-34.6 

6.1 

44.6 

18.0 

55173.39559 

-78.0 

5.7 

23.6 

16.2 

55197.34812 

-4.3 

4.8 

13.8 

13.3 

55221.51289 

-108.7 

4.3 

-3.2 

14.6 

55245.45360 

-216.0 

6.4 

-3.8 

14.8 

55358.13133 

-182.0 

5.4 

29.3 

13.3 

55512.49483 

56.7 

6.7 

46.8 

17.5 

55544.40921 

188.8 

6.8 

30.3 

17.6 



Table 3—Continued 


Epoch (MJD) 

RV (m s' 1 ) 

cr RV (m s x ) 

BS (m s x ) 

ubs (m s x ) 

55574.31534 

59.4 

5.4 

38.6 

13.4 

55577.32536 

64.5 

6.6 

18.5 

18.7 

55607.45310 

89.3 

6.3 

23.7 

15.0 

55637.15901 

37.6 

5.9 

65.4 

16.1 

55652.10131 

56.9 

4.9 

45.4 

11.6 

55691.21517 

-24.8 

6.0 

-0.7 

12.8 

55718.16090 

-89.4 

6.1 

0.4 

14.1 

55870.50090 

-141.0 

8.5 

1.5 

28.7 

55917.38294 

3.1 

5.2 

-17.7 

13.3 

55953.51961 

-42.6 

4.9 

33.6 

12.8 



Table 4. Relative radial velocities and bisector span of HD 216536 


Epoch (MJD) 

RV (m s' 1 ) 

0 -rv (m s x ) 

BS (m s x ) 

obs (m s *) 

53544.38601 

-43.1 

24.5 

-39.4 

17.4 

54365.34050 

24.9 

23.7 

64.5 

14.5 

54658.32385 

32.8 

25.1 

45.0 

30.8 

55370.39329 

53.7 

23.6 

-5.7 

15.4 

55470.35085 

-63.9 

23.8 

-12.0 

13.1 

55515.18324 

-44.2 

24.6 

-41.7 

12.4 

55518.18400 

-30.3 

24.0 

8.9 

14.9 

55728.36984 

8.8 

24.0 

26.5 

17.3 

55733.37812 

41.2 

23.7 

27.8 

13.5 

55740.37428 

4.0 

24.3 

-5.8 

17.8 

55753.32927 

-13.2 

23.8 

13.9 

14.2 

55758.31630 

-25.3 

23.9 

-4.8 

16.5 

55763.29701 

-21.3 

24.2 

6.2 

15.0 

55768.28240 

-54.0 

23.8 

20.3 

13.1 

55775.27116 

-53.2 

23.7 

19.6 

11.5 

55780.22843 

-68.7 

23.9 

5.6 

16.9 

55793.43019 

-57.6 

23.8 

34.3 

14.8 

55801.20164 

-25.2 

23.8 

-2.5 

13.3 

55810.16312 

-21.8 

23.8 

-15.7 

13.9 

55816.12715 

-40.3 

24.2 

1.0 

17.9 

55817.38332 

-8.9 

23.8 

-11.3 

14.9 

55822.36050 

4.5 

23.7 

27.5 

13.7 

55825.14133 

34.6 

23.9 

33.1 

14.6 

55826.14602 

49.4 

23.8 

29.2 

12.4 

55829.10961 

58.6 

23.7 

-13.9 

13.6 

55829.34425 

65.0 

23.7 

-7.1 

12.1 

55830.11034 

75.2 

23.7 

-3.2 

13.3 

55840.31604 

66.9 

23.8 

42.7 

11.3 

55849.31763 

5.4 

23.9 

25.4 

11.8 

55858.05332 

38.7 

23.9 

46.9 

12.3 

55888.21126 

-23.4 

23.7 

41.2 

12.1 



Tabic 4—Continued 


Epoch (MJD) 

RV (m s' 1 ) 

o-rv (m s x ) 

BS (m s x ) 

o-rs (m s x ) 

55903.14193 

-44.3 

24.0 

17.1 

17.8 

55924.08639 

-5.2 

23.9 

9.7 

14.1 

55932.06380 

-49.4 

23.7 

15.1 

10.8 

56073.43822 

-42.4 

23.8 

34.7 

17.5 

56080.43409 

-21.2 

23.7 

16.6 

12.4 

56179.39309 

25.8 

23.8 

23.8 

16.0 

56196.11293 

11.0 

23.7 

7.5 

13.4 

56210.06818 

-64.3 

23.7 

4.5 

11.7 

56228.25215 

-56.5 

23.9 

-7.7 

11.8 

56261.14362 

13.6 

23.6 

17.5 

11.6 

56262.16696 

34.9 

23.6 

12.1 

12.9 

56282.10282 

28.3 

23.9 

25.6 

15.9 

56287.10367 

4.6 

24.3 

40.7 

21.4 

56288.09887 

12.8 

23.8 

35.0 

11.6 

56291.08204 

44.2 

23.8 

17.2 

12.4 

56300.06651 

38.7 

23.8 

31.8 

13.2 
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Table 5. Orbital parameters of BD+49 828 b, HD 95127 b and HD 21653 b. 


Parameter 

BD+49 828 b 

HD 95127 b 

HD 21653 b 

p (DAYS) 

25901?“ 

482l(j 

148.61^ 

T 0 (MJD) 

554701?“ 

532001“ 

535871“ 

K (m s -1 ) 

18.8l®;g 

116lg 6 

501? 

e 

0-35lo;?o 

0 11+ 0 - 15 
u - ii -0.06 

0.381°;“ 

uj (deg) 

170+ 32 

1 ' u -30 

40+ 37 

^ u -40 

2701“ 

m 2 sin i (Mj) 

1 fi+0-4 
i - u - 0.2 

5 01 + 0 - 61 

J - ui -0.44 

1.471°;?° 

a (AU) 

A O+0.32 
^• Z - 0.2 

i 90 +O.OI 
^■• 20 - 0.01 

0.6091°;°“ 

V 0 (m s -1 ) 

1 1+ 1 ' 4 

—10.5111 

-4.91°;? 

Vid 

1.35 

1.14 

1.29 

o- RV (ms -1 ) 

11.6 

50.9 

23.0 

jitter (ms -1 ) 

4.441?;“ 

47.5 

17.9 

iVobs 

42 

41 

47 



Table 6. Basic parameters of planets around stars at various stages of stellar evolution 
(from exoplanets.eu). Mean and median values as well as dispersion are presented for every 

set of data. See text for details. 


Parameter M* range 

dwarfs 

subgiants 

giants 

bright giants 


mean a 

median 

mean a 

median 

mean a 

median 

mean a 

median 


N planet 

1-2 M 0 

267 


73 


23 


9 



all 

458 


77 


41 


14 



1-2 M 0 (RV) 

138 


60 


22 


9 



all(RV) 

236 


64 


40 


14 


log(g) 

all 

4.378 

0.009 

3.540 

0.037 

2.574 

0.039 

1.594 

0.077 



4.377 


3.394 


2.594 


1.663 



1-2 M e 

4.285 

0.008 

3.533 

0.037 

2.533 

0.054 

1.529 

0.100 



4.282 


3.386 


2.501 


1.521 



all(RV) 

4.392 

0.013 

3.485 

0.038 

2.564 

0.038 

1.594 

0.077 



4.392 


3.369 


2.582 


1.663 



1-2 M 0 (RV) 

4.296 

0.012 

3.473 

0.038 

2.514 

0.053 

1.529 

0.100 



4.300 


3.350 


2.501 


1.521 


M */M 0 

all 

1.039 

0.011 

1.428 

0.027 

1.870 

0.090 

1.464 

0.120 



1.031 


1.450 


1.900 


1.395 



1-2 M e 

1.179 

0.009 

1.453 

0.025 

1.499 

0.060 

1.457 

0.084 



1.140 


1.470 


1.500 


1.400 



all(RV) 

0.997 

0.015 

1.446 

0.031 

1.885 

0.091 

1.464 

0.120 



1.030 


1.475 


1.900 


1.395 



1-2 M e (RV) 

1.140 

0.010 

1.478 

0.028 

1.507 

0.062 

1.457 

0.084 



1.100 


1.480 


1.500 


1.400 


[Fe/H] 

all 

0.071 

0.010 

0.086 

0.022 

-0.032 

0.032 

-0.247 

0.058 



0.090 


0.120 


-0.052 


-0.255 



1-2 M 0 

0.114 

0.011 

0.099 

0.020 

-0.012 

0.045 

-0.222 

0.055 



0.140 


0.140 


-0.030 


-0.260 



all(RV) 

0.102 

0.014 

0.093 

0.026 

-0.036 

0.033 

-0.247 

0.058 



Table 6—Continued 


Parameter M* range 

dwarfs 

subgiants 

giants 

bright giants 


mean a 

median 

mean a 

median 

mean a 

median 

mean a 

median 




0.140 


0.140 


-0.056 


-0.255 



1-2 M 0 (RV) 

0.178 

0.210 

0.012 

0.109 

0.140 

0.023 

-0.018 

-0.060 

0.046 

-0.222 

-0.260 

0.055 

m psini / m j 

all 

2.046 

0.955 

0.147 

2.702 

1.800 

0.399 

6.391 

4.500 

0.987 

7.798 

6.300 

1.186 


1-2 M 0 

2.459 

1.308 

0.212 

2.705 

1.800 

0.416 

4.557 

3.200 

0.592 

7.688 

7.800 

0.854 


all(RV) 

2.250 

1.135 

0.200 

2.956 

1.850 

0.467 

6.533 

4.900 

1.001 

7.798 

6.300 

1.186 


1-2 M 0 (RV) 

2.532 

1.658 

0.256 

2.977 

1.850 

0.492 

4.731 

3.260 

0.592 

7.688 

7.800 

0.854 

e 

all 

0.162 

0.072 

0.010 

0.185 

0.160 

0.018 

0.187 

0.110 

0.031 

0.220 

0.167 

0.048 


1-2 M 0 

0.165 

0.070 

0.013 

0.178 

0.157 

0.018 

0.177 

0.090 

0.045 

0.207 

0.140 

0.068 


all(RV) 

0.244 

0.201 

0.014 

0.212 

0.185 

0.019 

0.190 

0.111 

0.032 

0.220 

0.167 

0.048 


1-2 M 0 (RV) 

0.260 

0.225 

0.018 

0.205 

0.179 

0.019 

0.182 

0.095 

0.047 

0.207 

0.140 

0.068 

a/AU 

all 

1.116 

0.094 

0.193 

1.331 

1.310 

0.111 

1.476 

1.200 

0.172 

1.487 

1.470 

0.085 


1-2 M 0 

1.442 

0.095 

0.324 

1.358 

1.310 

0.115 

1.119 

1.160 

0.113 

1.517 

1.540 

0.119 


all(RV) 

1.506 

0.983 

0.108 

1.568 

1.450 

0.112 

1.511 

1.235 

0.173 

1.487 

1.470 

0.085 
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Table 6—Continued 


Parameter M* range 

dwarfs 

subgiants 

giants 

bright giants 


mean 

median 

a mean 

median 

a mean 

median 

cr mean cr 

median 

1-2 M 0 (RV) 

1.725 

1.188 

0.149 1.617 

1.515 

0.115 1.166 

1.180 

0.107 1.517 0.118 

1.540 


Table 7. Parsons correlation coefficients between stellar hosts mass and total mass of their 
planetary systems. Upper values: all systems for which at least masses are available; lower 
values: systems with all data listed in Table 6 available. All correlation coefficients are 

smaller than respective critical values. 


group 


M 

PS vs. M* 




N systems 



1-2M© 

all 

1-2M© (RV) 

all (RV) 

1-2M© 

all 

1-2M© (RV) 

all (RV) 

dwarfs 

0.222 

0.132 

0.194 

0.234 

256 

432 

100 

172 

0.137 

0.172 

0.202 

0.223 

220 

367 

99 

167 

subgiants 

0.128 

0.133 

0.110 

0.109 

68 

73 

55 

59 

0.119 

0.113 

0.095 

0.099 

64 

68 

53 

57 

giants 

0.101 

0.409 

0.057 

0.400 

21 

40 

20 

39 

0.182 

0.480 

0.144 

0.469 

20 

37 

19 

36 

bright giants 

0.114 

0.604 

0.114 

0.604 

9 

16 

9 

16 

0.114 

0.606 

0.114 

0.606 

9 

14 

9 

14 

all 

0.204 

0.234 

0.164 

0.402 

354 

56#| 

184 

286 

0.181 

0.380 

0.172 

0.407 

313 

486 

180 

274 
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Fig. 1.— The PTPS sample. Presented are 838 stars with detailed spectroscopic parameters 
available (Zielinski et al. 2012 Niedzielski cl al. in prep., Deka et ah in prep.). The three 
stars presented in this paper are marked as black rectangles. Color coded is the amplitude 
of p-mode oscillations estimated from the scaling relations of Kjcldsen & Bedding; (1995). 
Evolutionary tracks of Bertelli et ah (2008) are presented for three metallicities: [Fe/H]=0.0 
(solid line), -0.30 (dashed line) and 0.3 (dotted line) for illustration purposes. 
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Fig. 2.— Radial velocities of BD+49 828 as in Table [2] plotted as a function of time together 
with the solution detailed in Table |5l 
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Fig. 3.— Radial velocities of HD 95127 as in Table [4] plotted as a function of time together 
with the solution detailed in Table |5l 
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Fig. 4.— Radial velocities of HD 216536 as in Table [4] plotted as a function of time together 
with the solution detailed in Table |5l 
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Phase 


Fig. 5.— Radial velocities of BD+49 828 plotted as a function of orbital phase together with 
the solution detailed in Table [5} 
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Fig. 6.— Radial velocities of HD 95127 plotted as a function of orbital phase together with 
the solution detailed in Table [5] 
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Phase 


Fig. 7.— Radial velocities of HD 216536 plotted as a function of orbital phase together with 
the solution detailed in Table [5} 
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Fig. 


Lomb-Scargle periodograms for BD+49 828. From top to bottom: HET RV data, 


BS, WASP photometry (Butters et al. 2010), /h«, Ipe, and RV post-fit residua. The low- 
signal peak in the periodogram for WASP photometry is most certainly an artifact as the 
available data cover only a fraction of the orbital period present in the RV data - see Sect. 
4.1 for discussion. 
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Fig. 9.— Lomb-Scargle periodograms for HD 95127. From top to bottom: HET RV data, 


BS, Hipparcos photometry (ESA 1997), Jh«, /pe, and RV post-fit residua. 
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18 



Period [days] 


Fig. 10.— 
BS, NSVS 


Lomb-Scargle periodograms for HD 216536. From top to bottom: HET RV data, 
photometry (Wozniak et al. 2004), Iu a , Ip e , and RV post-fit residua. 
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Fig. 11.— Orbital separation - log(a/R*) for planets orbiting hosts at various evolutionary 
stages - log(g), all data from exoplanets.eu. Grey symbols - all planets, black symbols - 
planets around 1 — 2 M 0 stars, red symbols - new planets presented here. Red and blue 
lines - minimum distance to avoid engulfment for a 1.5 M 0 star and 1 Mj planet system and 
minimum orbit above which a planet is not sensitive to tidal interaction with the stellar host 
from Villaver et ah (2014). See text for discussion. 








